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When quarks and gluons tend to form a dense medium, like in high energy or/and heavy-ion 
collisions, it is interesting to ask the question which are the relevant degrees of freedom that Quantum 
Chromodynamics predict. The present notes correspond to two lectures given at Zakopane in the 
(rainy) summer of 2006, where this question is adressed concretely in two cases, one in the QCD 
regime of weak coupling, the other one at strong coupling. Each case corresponds to the study of an 
elusive but dynamically important transient phase of quarks and gluons expected to appear from 
Quantum Chromodynamics during high energy collisions. 

In lecture I, we examine the dynamical phase space of gluon transverse momenta near the so-called 
"saturation" phase including its fluctuation pattern. "Saturation" is expected to appear when the 
density of gluons emitted during the collision reaches the limit when recombination effects cannot be 
neglected, even in the perturbative QCD regime. We demonstrate that the gluon-momenta exhibit 
a nontrivial clustering structure, analoguous to "hot spots", whose distributions are derived using 
an interesting matching with the thermodynamics of directed polymers on a tree with disorder and 
its "spin-glass" phase. 

In lecture II, we turn towards the non-perturbative regime of QCD, which is supposed to be 
relevant for the description of the transient phase as quark-gluon plasma formed during heavy-ion 
collisions at very high energies. Since there is not yet an available field-theoretical scheme for non 
perturbative QCD in those conditions, we study the dynamics of strongly interacting gauge-theory 
matter (modelling quark-gluon plasma) using the AdS/CFT duality between gauge field theory at 
strong coupling and a gravitational background in Anti-de Sitter space. The relevant gauge theory 
is a-priori equipped with TV = 4 supersymmetries, but qualitative results may give lessons on this 
issue. As an explicit example, we show that perfect fluid hydrodynamics emerges at large times as 
the unique nonsingular asymptotic solution of the nonlinear Einstein equations in the bulk. The 
gravity dual can be interpreted as a black hole moving off in the fifth dimension. 
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LECTURE I: QCD NEAR SATURATION: A SPIN-GLASS STRUCTURE 



Introduction 



1. Saturation in QCD is expected to occur when parton densities inside an hadronic target are so high that their wave- 
functions overlap. This is expected from the rapidity Y — log(W^) evolution of deep-inelastic scattering amplitudes 
governed by the Balitsky Fadin Kuraev Lipatov (BFKL) kernel lj. The BFKL evolution equation is such that the 
number of gluons of fixed size increases exponentially and would lead without modification to a violation of unitarity. 
By contrast, the renormalisation-group evolution equations following Dokshitzcr, Gribov and Lipatov, Altarelli and 
Parisi (DGLAP) Q explains the evolution at fixed Y as a function of the hard scale Q 2 . It leads to a dilute system of 
asymptotically free partons. As schematized in Fig.l, the transition to saturation 0,0 is characterized by a typical 
transverse momentum scale Q S (Y), depending on the overal rapidity of the reaction, when the unitarity bound is 
reached by the BFKL evolution of the amplitude. The two-dimensional plot showing the two QCD evolution schemes 
and the transition boundary to saturation are represented in Fig.l. 



c-4 
< 



M 
O 
-J 



$5 



CGC? 




BFKL 



DGLAP 
> 



L = Log(Q*2) 



FIG. 1: Schematic view of the transition region to saturation. The DGLAP and BFKL evolution ranges are displayed, together 
with the saturation region where the density bounds are reached. 



The problem we address here is the characterization of the gluon-momentum distribution near saturation. Our 
aim is to understand the transverse-momenta spectrum of the gluons which are generated by the BFKL evolution 
in rapidity, i.e. as characterising a transient QCD phase structure near saturation as a whole. The new material 
contained in this lecture comes from Ref. ||. 

As a guide for the further developments, the basic structure underlying the transition to saturation can be un- 
derstood in terms of traveling waves. If at first one neglects the fluctuations (in the mean-field approximation), the 
effect of saturation on a dipole-target amplitude is described by the nonlinear Balitsky-Kovchegov p (BK) equation, 
where a nonlinear damping term adds to the BFKL equation. As shown in Q, this equation falls into the universality 
class of the Fisher and Kolmogorov Petrovsky Piscounov (F-KPP) nonlinear equation Q which admits asymptotic 
traveling-wave solutions. The exponential behaviour of the BFKL evolution quickly enhances the effects of the tail 
towards a region where finally the nonlinear damping regulates both the traveling-wave propagation and structure. 

Indeed, one of the major recent challenges in QCD saturation is the problem of taking into account the role of 
fluctuations, i.e. the structure of gluon momenta beyond the average. In these conditions it was realized for traveling 
waves and thus in the QCD case that the fluctuations may have a surprisingly large effect on the overall 
solution of the nonlinear equations of saturation. Indeed, a fluctuation in the dilute regime may grow exponentially 
and thus modify its contribution to the overall amplitude. Hence, in order to enlarge our understanding of the QCD 
evolution with rapidity, it seems important to give a quantitative description of the pattern of momenta generated by 
the BFKL evolution equations for the set of cascading dipoles (or, equivalently gluons) near saturation, which is the 
subject of the lecture. 



Technically speaking, we shall work in the leading order in 1/N C , where the QCD dipole framework is valid 
Moreover we will use the diffusive approximation of the 1-dimensional BFKL kernel. In fact, the phase structure 
appears quite rich already within this approximation scheme. Many aspects we will obtain show "universality" 
features and thus are expected to be valid beyond the approximations. 



Rapidity evolution of cascading QCD dipoles 

2. Let us start by briefly describing the QCD evolution of the dipole distributions [HHl^. Ilflj. 
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FIG. 2: BFKL cascading and saturation. The QCD branching process in the BFKL regime and beyond is represented along 
the rapidity axis (upper part). Its 2-dimensional counterpart in transverse position space is displayed at two different rapidities 
(lower part). The interaction region is represented by a shaded disk of size 1/Q. At rapidity Y\, the interaction still probes 
individual dipoles (or gluons), which corresponds to the exponential BFKL regime. There exists a smooth transition to a regime 
where the interaction only probes groups of dipoles or gluons, e.g. at rapidity Yg. This gives a description of the near- saturation 
region corresponding to a mean-field approximation, where correlations can be neglected. Further in rapidity, Y > Y2, other 
dynamical effects, such that merging and correlations appear. 



The structure of BFKL cascading describes a 2-dimensional tree structure of dipoles in transverse position space 
evolving with rapidity. Let us, for instance, focus on the rapidity evolution starting from one massive qq pair or onium 
see Fig|21 At each branching vertex, the wave function of the onium-projectile is described by a collection of 
color dipoles. The dipoles split with a probability per unit of rapidity defined by the BFKL kernel [l| 

/C(v, w; z) = a s N c /n (v ~ w) (1) 
(v — z) (z — w) 

describing the dissociation vertex of one dipole (v, w) into two dipoles at (v, z) and (z, w), where v, w, z are arbitrary 
2-dimensional transverse space coordinates. 

As an approximation of the 2-dimensional formulation of the BFKL kernel JU obtained when one neglects the 
impact-parameter dependence, we shall restrict our analysis in the present paper to the 1-dimensional reduction of 
the problem to the transverse-momenta moduli hi of the cascading gluons. After Fourier transforming to transverse- 
momentum space, the leading-order BFKL kernel defining the rapidity evolution in the 1-dimensional approxima- 
tion is known 0] to act in transverse momentum space as a differential operator of infinite order 



x (_$) = 2^(l)-^(-$)-V(l + ft) 



(2) 



where I = logfc 2 and Y is the rapidity in units of the fixed coupling constant a s N c /ir. 

In the sequel, we shall restrict further our analysis to the diffusive approximation of the BFKL kernel. We thus 
expand the BFKL kernel to second order around some value 7 C 

X(7) ~ Xc + x'Al - 7c) + |Xc(7 ~ 7c) 2 = Au - A l7 + A 2l 2 , (3) 

where j c will be defined in such a way to be relevant for the near-saturation region of the BFKL regime. 

Within this diffusive approximation, it is easy to realize that the first term (Aq) is responsible for the exponential 
increase of the BFKL regime while the third term (A 2 ) is a typical diffusion term. The second term (A\ ) is a "shift" 
term since it amounts to a rapidity-dependent redefinition of the kinematic variables, as we shall see. 

In Eq.®, 7c is chosen in order to ensure the validity of the kernel © in the transition region from the BFKL 
regime towards saturation. Indeed, the derivation of asymptotic solutions of the BK equation Q leads to consider 
the condition 

X(7c) = 7c x'(7c) (4) 

whose solution determines 7 C . 

This condition applied to the kernel formula gives j c = \/ Aq/A 2 ~ 0.6275... and {Aq, A\, A 2 } « 
{9.55,25.56,24.26}. These numbers may appear anecdotic, but they fully characterize the critical parameters of 
the saturation transition, as we will realize later on. For different kernels, e.g. including next-leading log effects [l4| . 
they could be different, of course. But, then the traveling-wave solution will be in a different "universality class" in 
mathematical terms. 



Mapping to thermodynamics of directed polymers 

3. From the properties in transverse-momentum space and within the I-dimensional diffusive approximation J3jl, 
we already noticed that the BFKL kernel models boils down to a branching, shift and diffusion operator acting in 
the gluon transverse-momentum-squared space. Hence the cascade of gluons can be put in correspondence with a 
continuous branching, velocity-shift and diffusion probabilistic process, see Fig|2 whose probability by unit of rapidity 
is defined by the coefficients A; of 
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FIG. 3: Branching diffusion model for polymers. The coordinates xi(t) ■ ■ ■ xj{t) correspond to the random paths along the tree 
in the x,t phase-space. The oblique axis is for L = f3x + Ai/Ao t, which takes into account the "time-drift" in the mapping to 
the QCD problem. 

Let us first introduce the notion of gluon- momenta "histories" fej(y). They register the evolution of the gluon- 
momenta starting from the unique initial gluon momentum k(0) and terminating with the specific z-th momentum 
hi, after successive branchings. They define a random function of the running rapidity y, with < y < Y, the final 



rapidity range when the evolution ends up (say, for a given total energy). It is obvious that two different histories 
ki{y) and kj(y) are equal before the rapidity when they branch away from their common ancestor. 

In order to fromulate precisely the mapping to the polymer problem, we then introduce random paths Xi(t) using 
formal space and time coordinates x, t , < t < At which we relate to gluon-momenta "histories" as follows: 



t 

A~ 



log A;? (y) = -(3 {x t {t) - x(0)) + (A - A^y 



(5) 



where (A$ — Ai)y is a conveniently chosen and deterministic "drift term", x(0) is an arbitrarily fixed origin of an 
unique initial gluon and thus the same for all subsequent random paths. The random paths Xi(t) are generated by a 
continuous branching and Brownian diffusion process in space-time (cf. Fig|3J|. 

As we shall determine later on, the important parameter f3, which plays the role of an inverse of the temperature 
T for the Brownian movements of the polymer process, is given by 



1 

T 



13 = ^2A 2 /A 



(G) 



In fact, the relation JHl will be required by the condition that the stochastic process of random paths describes the 
BFKL regime of QCD near saturation. Another choice of [3 would eventually describe the same branching process 
but in other conditions. Hence the condition (jSJ) will be crucial to determine the QCD phase at saturation (within 
the diffusive approximation). 

Let now introduce the tree-by-tree random function defined as the partition function of the random paths system 



Z(t) = S^ e ~^ Xi( ^ = e ~P Xo+AlV 

i=l 



*f(y) cc 



x P(y) , 



(7) 



where ~52iLi ^iiv) = k 2 {y) is the event-by-event average over gluon momenta at rapidity y. Note that one has to 
distinguish i.e. the average made over only one event from ( • • • ), which denotes the average over samples (or 

events) . 

Z(t) is an event-by-event random function. The physical properties are obtained by averaging various observables 
over the events. Note that the distinction between averaging over one event and the sample-to-sample averaging 
appears naturally in the statistical physics problem in terms of "quenched" disorder: the time scale associated with 
the averaging over one random tree structure is much shorter than the one corresponding to the averaging over random 
trees. 

With these definitions, Z appears to be nothing else than the partition function for the model of directed polymers 
on a random tree [lfij ]. 

Let us now justify the connection of the directed-polymer properties with the description of the gluon-momentum 
phase near saturation by rederiving the known saturation features from the statistical model point-of-view. Using the 
known properties of the partition function of the polymer problem, one finds 



logQ2 = {logfc 2 



logZ)-A 1 Y 



(8) 



which in fact matches exactly the asymptotic expansion found in for the saturation scale. 

In the same way, the solution of the statistical-physics problem allows to derive the event-by-event spectrum of the 
free energy log Z(t) of the system around its average. From 0, one gets 



Af(k 2 ,Y) ~V(\ogZ- (logZ)) oc log 



r>2 



exp 




(9) 



which is the well-known geometrical scaling property, empirically found in Ref.|l7| and theoretically derived in Q 
from the BK equation for the dipole amplitude J\f(k 2 , Y). 

Both properties ©,© prove the consistency of the model with the properties expected from saturation. We shall 
then look for other properties of the cascading gluon model. It is important to realize that this consistency fails for 
a different choice of the parameter j3 different from ©. This justifies a-posteriori the identification of the equivalent 
temperature of the system in the gluon/polymer mapping framework. 



The Spin-Glass Phase of gluons 



4. Let us now come to the main new results concerning the determination of structure of the gluon-momentum phase 
at saturation. 

The striking property of the directed polymer problem on a random tree is the spin-glass structure of the low 
temperature phase. As we shall see this will translate directly into a specific clustering structure of gluon transverse 
momenta in their phase near the "unitarity limit", see Fig0J 
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FIG. 4: The clustering structure 0} gluons near saturation. The drawing represents the si • • • S4 clusters near momenta k s i ■ ■ ■ k S 4. 
They branch either near t -C 1 or (At — t) -C 1, where At is the total amount of time evolution. 

Following the relation (j^J , one is led to consider the polymer system at temperature T with 

T c -T B c 

-V-l-^=l-7c (10) 

where 7 C is the critical exponent defined by (J3J. We are thus naturally led to consider the low-temperature phase 
(T < T c ), at some distance T c — T from the critical temperature T c =i/V2. In the language of traveling waves [lfij . 
this corresponds to a "pulled- front" condition with "frozen" and "universal" velocity and front profile. 

As derived in [l5j |. the phase space of the polymer problem is structured in "valleys" which are in the same 
universality class as those of the Random Energy Model (REM) and of the infinite range Sherrington-Kirkpatrick 
(SK) model [H- 

Translating these results in terms of gluon-momenta moduli, the phase space landscape consists in event- by- event 
distribution of clusters of momenta around some values h 2 si = \j{rii) ^2 i£si kf, where rii is the cluster multiplicity. 
The probability weights to find a cluster Sj after the whole evolution range Y is defined by 

W sl = ^g^, (11) 

where the summation in the numerator is over the momenta of gluons within the sf l cluster (see Fig0J. The normalized 
distribution of weights W S i thus allows one to study the probability distribution of clusters. The clustering tree 
structure, (called "ultrametric" in statistical mechanics) is the most prominent feature of spin-glass systems |2l|. 

Note again that, for the QCD problem, this property is proved for momenta in the region of the "unitarity limit", or 
more concretely in the momentum region around the saturation scale. This means that the cluster average-momentum 
is also such that k 2 si = (D(Q 2 S ). Hence the clustering structure is expected to appear in the range which belongs to the 
traveling-wave front [?J or, equivalently, of clustering with finite fluctuations around the saturation scale. 
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FIG. 5: The probability distribution of overlaps Tl(y). The figure (the simulation is by courtesy from [Tsl| . using the method 
of Ref.|23() is drawn both for the theoretical QCD value 7 C = 0.7, and for 7 C = 0.3 for comparison. The statistics used for the 
simulation is 10 s events in 10 3 bins for y c = 0.3 and 3.25 10 6 events in 10 2 bins for 7 C = 0.7. 

In order to quantify the cluster structure, one may introduce a well-known "overlap function" in statistical physics 
of spin-glasses |2lJ. Translating the definitions (cf. 15]) in terms of the QCD problem, one introduces an event- 
by-event indicator of the strength of clustering which is built in from the weights (fill) , namely y = J2 S i ^ nc 
non-trivial probability distribution of overlaps H(y) possesses some universality features, since it is identical to the 
one of the REM and SK models and shares many qualitative similarities with other systems possessing a spin-glass 
phase USUI]- Exemples are given in FigEl 

The rather involved probability distribution H(y) is quite intringuing. It possesses a priori an infinite number of 
singularities at y — 1/n, n integer. It can be seen when the temperature is significantly lower from the critical value, 
see for instance the curve for 1 — T/T c = .7 in Fig|3] However, the predicted curve for the QCD value 1 — T/T c ~ .3 
is smoother and shows only a final cusp at W s = 1 within the considered statistics. It thus seems that configurations 
with only one cluster can be more prominent than the otherwise smooth generic landscape. However, it is also a 
"fuzzy" landscape since many clusters of various sizes seem to coexist in general. 

Summary of lecture I 

5. We investigated the landscape of transverse momenta in gluon cascading around the saturation scale at asymptotic 
rapidity. Limiting our study to a diffusive 1-dimensional modelization of the BFKL regime of gluon cascading, we 
make use of a mapping on a statistical physics model for directed polymers propagating along random tree structures 
at fixed temperature. We then focus our study on the region near the unitarity limit where information can be 
obtained on saturation, at least in the mean-field approximation. Our main result is to find a low-temperature spin- 
glass structure of phase space, characterized by event-by-event clustering of gluon transverse momenta (in modulus) 
in the vicinity of the rapidity-dependent saturation scale. The weight distribution of clusters and the probability of 
momenta overlap during the rapidity evolution are derived. 

Interestingly enough the clusters at asymptotic rapidity are branching either near the beginning ("overlap 0" or 
y/Y <C 1) or near the end ("overlap 1" or 1 — y/Y <C 1) of the cascading event. The probability distribution of 
overlaps is derived and shows a rich singularity structure. 

On a phenomenological ground, it is remarkable that saturation density effects are not equally spread out on the 
event-by-event set of gluons; our study suggests that there exists random spots of higher density whose distribution 
may possess some universality properties. In fact it is natural to expect this clustering property to be present not only 
in momentum modulus (as we could demonstrate) but also in momentum azimuth-angle. This is reminiscent of the 
"hot spots" which were some time ago Q advocated from the production of forward jets in deep-inelastic scattering 
at high energy (small-x). The observability of the cluster distribution through the properties of "hot spots" is an 
interesting possibility. 




LECTURE II: A PERFECT FLUID FROM GAUGE/GRAVITY DUALITY 

Introduction 



6. From the first years of the running of heavy-ion collisions at RHIC, evidence has been advocated that various 
observables are in good agreement with models based on hydrodynamics [2j| and with quark-gluon plasma (QGP) in 
a strongly coupled regime [2||. To a large extent it seems that the QGP behaves approximately as a perfect fluid as 
was first considered in (S^. A schematic view of the theoretical expectations in given in FigO It is a challenge in QCD 
to derive from first principles the properties of the dynamics of a strongly interacting plasma formed in heavy-ion 
collisions and in particular to understand why the perfect-fluid hydrodynamic equations appear to be relevant. 



FIG. 6: Scenario for the quark-gluon plasma( QGP) formation. After a pre-equilibrium stage in a heavy-ion collision, probably 
governed by a weak-coupling but dense regime, the QGP is formed with local equilibrium and hydrodynamic properties. 

Even if the experimental situation is still developing and rather complex, it is worth simplifying theproblem in order 
to be able to attack it with appropriate theoretical tools. Recently the AdS/CFT correspondence j23,|2i| emerged as 
a new approach to study strongly coupled gauge theories. This has been largely worked out in the supersymmetric 
case and in particular for the conformal case of J\f = 4 super Yang-Mills theory (SYM). Interestingly enough, since 
the QGP is a deconfined and strongly interacting phase of QCD we could expect that results for the nonconfining 
M = 4 theory may be relevant or at least informative on the unknown strong coupling QCD problem. We will start 
from this assumption in our work. 

In this lecture we focus on the spacetime evolution of the gauge theory (4d) energy- momentum tensor, and derive 
its asymptotic behaviour from the solutions of the nonlinear Einstein equations of the gravity dual. 

Imposing the absence of curvature singularities in the gravity dual, we will show that, in the boost invariant setting 
(as in (13), perfect fluid hydro dynamics emerges from the AdS/CFT solution at large times. The new material 
contained here comes from Refs. I3(M. 



7. As an introduction to our lecture, let us briefly recall some aspects of the String/Gauge Duality. The AdS/CFT 
correspondence [i^ has many interesting formal and physical facets. Concerning the aspects which are of interest for 
our problem, it allows one to find relations between gauge field theories at strong coupling and string gravity at weak 
coupling in the limit of large number of colours (N c — >oo). It can be examined quite precisely in the AdSs/CFT4 case 
which conformal field theory corresponds to SU(N) gauge theory with J\f = 4 supersymmetries. 

Let us recall the canonical derivation leading to the AdSs background , see Fig0 One starts from the (super)gravity 
classical solution of a system of N Z?3-branes in a 10 — D space of the (type IIB) superstrings. The metrics solution 
of the (super)Einstein equations read 
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where the first four coordinates are on the brane and r corresponds to the coordinate along the normal to the branes. 
In formula (|12|) . one defines 



f = ! + ?-■ R = ^g 2 YM a' 2 N , 



(13) 



where g YM N is the 't Hooft- Yang-Mills coupling and a' the string tension. One considers the limiting behaviour 
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FIG. 7: Schematic view of the Gauge/String Duality. Left: The string background is the Anti-de Sitter space (AdS); Right: 
the gauge theory is a Conformal field theory (CFT) on the 4-d N-branes. 

considered by Maldacena, where one zooms on the neighbourhood of the branes while in the same time going to the 
limit of weak string slope a' . The near-by space-time is thus distorted due to the (super) gravitational field of the 
branes. One goes to the limit where 



R fixed 



r(- 0) 

This, from the second equation of Ijl3|) obviously implies 

a' , g 2 YM N ~ i- 



z fixed 



(14) 



(15) 



i.e. both a weak coupling limit for the string theory and a strong coupling limit for the dual gauge field theory. By 
reorganizing the two parts of the metrics one obtains 



ds 2 = -^(-dt 2 + dx 2 + dz 2 ) + R 2 dn 5 



(16) 



1-3 



which corresponds to the AdS?, x S§ background structure, 55 being the 5-sphere. More detailed analysis shows 
that the isometry group of the 5-sphere is the geometrical dual of the Af = 4 supersymmetries. More intricate is the 
quantum number dual to N c , the number of colours, which is the invariant charge carried by the Ramond-Ramond 
form field. 



Bjorken hydrodynamics 



8. Coming back to the physical world, a model of the central rapidity region of heavy-ion reactions based on hydro- 
dynamics was pioneered in [22| and involved the assumption of boost invariance. Our goal is to study the dynamics 
of strongly interacting gauge-theory matter assuming boost invariance. 

We will be interested in the spacetime evolution of the energy-momentum tensor T M „ of the gauge-theory matter. 
It is convenient to introduce proper-time (r) and space-rapidity (y) coordinates in the longitudinal position plane: 

x° = t coshy x 1 = t siiiliy . (17) 

In these coordinates, all components of the energy momentum tensor can be expressed (see |30j ]') in terms of a single 
function /(r): 



T 



( f(r) \ 

-^/(t)-t 2 /(t) 

f(T)+%T-fcf(T) 

V f(T) + §T&f(T)J 



(18) 



where the matrix T M „ is expressed in (r, y, x\,X2) coordinates. 
Furthermore the function fir) is constrained to verify 

/(r)>0 /'(r)<0 r/'(r) > -4/(r) . (19) 

The dynamics of the gauge theory should pick a specific /(t). A perfect fluid or a fluid with nonzero viscosity and/or 
other transport coefficients will lead to different choices of /(t). 

We thus address the problem of determination of the function /(r) from the AdS/CFT correspondence. Let us first 
describe two distinct cases of physical interest: 

For a perfect fluid (Bjorken hydrodynamics) /(r) ~ l/i~3, while for a "free streaming case" expected just at the 
beginning of the interaction , /(t) ~ 1 jr. In the following we introduce a family of /(r) with the large r behaviour 
of the form 

/(r)~r- s . (20) 

Boost-invariant geometries 

9. The most general form of the bulk metric respecting boost-invariance can be written 

_ e a{r,z) dT 2 +T 2 e b{r,z) d 2 + e c{r,z) dx 2 dz 2 

ds 2 y - ^ + — . (21) 

The three coefficient functions can be (non trivially) derived from the Einstein equations 

R tiV -~g^R-&g^ = Q , (22) 

in the asymptotic limit where r — > oo. Interestingly enough, they depend only on the scaling variable v = z/t s / 4 , 
where s labels the one-parameter family of solutions 120(1 . 
After quite painful calculations, the solution reads |3(| 

a(v) = A(v) - 2m(v) b{v) = A(v) + (2s - 2)m(u) c(v) = A{v) + (2 - s)m{v) (23) 

where 

A(v) = l - (log(l + AGO v 4 ) + log(l - A( s ) v 4 )) m(v) = (log(l + A(s) v 4 ) - log(l - A(s) v 4 )) (24) 

with 




FIG. 8: The curvature scalar d\ 2 . JH 2 is calculated as a function of w = n/A(s) 1//4 for the perfect fluid case s = 4/3 (solid line), 
s = 4/3 - 0.1 (dotted line) and s = 4/3 + 0.2 (dashed line). 



Specializing first to the perfect fluid case, this gives rise to the following asymptotic geometry 

( l — e ° z A \ 

z2 ds " = - \ !/Z dT ' 2 + f 1 + f 7^) W + + dz ' W 
1 + 3 7*73 

Remarkably enough this geometry can be identified (in suitable metrics) to be a moving Black Hole, which evolves in 
the fifth dimension z. 

For the free streaming case, one finds 

9,9 , V 4 . 1-2 vg , V 4 . 1+2 , , , V 4 . 1 , f 4 x 1 9 , 9 , « 4 s ifa/g V 4 s W2 „ „ 

ds 2 = -(i + _ _)-b-(i -)—2—dt 2 + 1 + — 2 1 =)*T 2 dy 2 + (i + _)-3-(i — dx\ +dz 2 , 27 

y/8 V8 V8 VS V8 V8 

which is qualitatively different from the perfect fluid case, in particular it displays singularities or zeroes at v 4 = \/8 
in all coefficients. More generally, it is possible to show that the perfect fluid case is the only one free of physical 
singularities, namely singularities which cannot be removed by a change of coordinates. In order to check this feature, 
we considered the metric-invariant curvature scalar 

m 2 - R^R^p . (28) 

As an illustration, we represent this property in FigEI where the value of %\ 2 is studied as a function of the distance 
from the horizon, for s values at the perfect fluid point and very near-by values. Let us add some comments on the 
specific features of our approach and results. We concentrate on looking for solutions of the full nonlinear Einstein 
equations. It would be interesting to confront this approach with the linearization methods of refs. |3ll |. In particular 
viscosity terms are expected to appear in the study of subasymptotic terms. Note that the possibility of black hole 
formation in the dual geometry has been argued in ref. |32l. More specifically, the geometry of a brane moving w.r.t. 
a black hole background has been advocated in ref. [33j for the dual description of the cooling and expansion of a 
quark-gluon plasma. In our case we could interpret the solution (|26|l as a kind of 'mirror' situation in terms of a black 
hole moving off from the AdS boundary. 



Summary of Lecture II 



10. We have introduced a general framework for studying the dynamics of matter (plasma) in strongly coupled gauge 
theory using the AdS/CFT correspondence for the J\f = 4 SYM theory. We constructed dual geometries for given 
4-dimensional gauge theory energy-momentum tensor profiles. Further imposing boost-invariant dynamics inspired by 
the Bjorkcn hydrodynamic picture, we have found the corresponding asymptotic solutions of the nonlinear Einstein 
equations. Among the family of asymptotic solutions, the only one with bounded curvature scalars is the gravity 
dual of a perfect fluid through its energy-momentum tensor profile. This selected nonsingular solution, given by the 
metric (|26[) . corresponds to a black hole moving off in the 5th dimension as a function of the physical proper time. 
As an application of this framework, we can obtain the thermalization time of the perfect fluid, which describes 
the decay back to equilibrium of a scalar excitation of the perfect fluid out of equilibrium, by computation of the 
quasi-normal modes of the moving Black Hole. In some sense, the moving Black Hole is a quite stable geometric 
configuration. We conjecture that it may represent, through the Gauge/Gravity duality, a powerful "attractor" for 
the QGP evolution, or even perhaps for more general evolution of a strongly coupled system of quarks and gluons. 
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